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In brief 
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individuals vaccinated against SARS- 
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of cross-reactive antibodies to seasonal 
coronaviruses and also maps the regions 
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SUMMARY 


In this study we profiled vaccine-induced polyclonal antibodies as well as plasmablast-derived mAbs from 
individuals who received SARS-CoV-2 spike mRNA vaccine. Polyclonal antibody responses in vaccinees 
were robust and comparable to or exceeded those seen after natural infection. However, the ratio of binding 
to neutralizing antibodies after vaccination was greater than that after natural infection and, at the mono- 
clonal level, we found that the majority of vaccine-induced antibodies did not have neutralizing activity. 
We also found a co-dominance of mAbs targeting the NTD and RBD of SARS-CoV-2 spike and an original 
antigenic-sin like backboost to spikes of seasonal human coronaviruses OC43 and HKU1. Neutralizing activ- 
ity of NTD mAbs but not RBD mAbs against a clinical viral isolate carrying E484K as well as extensive changes 
in the NTD was abolished, suggesting that a proportion of vaccine-induced RBD binding antibodies may pro- 


vide substantial protection against viral variants carrying single E484K RBD mutations. 


INTRODUCTION 


Understanding of the innate and adaptive immune responses to 
severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) 
has progressed rapidly since the beginning of the coronavirus 
disease 2019 (COVID-19) pandemic (Carvalho et al., 2021). Poly- 
clonal antibody responses against the spike protein of the virus 
in serum, and to a lesser degree also at mucosal surfaces, 
have been well characterized with respect to their kinetics, bind- 
ing capacity, and functionality (Grandjean et al., 2020; Isho et al., 
2020; lyer et al., 2020; Ripperger et al., 2020; Seow et al., 2020; 
Wajnberg et al., 2020). Similarly, encouraging data have been 
published on both the plasmablast response and the memory 
B cell response induced by SARS-CoV-2 infection (Dan et al., 
2021; Gaebler et al., 2020; Guthmiller et al., 2021; Huang et al., 
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2021; Robbiani et al., 2020; Rodda et al., 2021; Wilson et al., 
2020). The immune responses to SARS-CoV-2 vaccination, 
including to mRNA-based vaccines, are less well studied since 
these vaccines only became available in the last months of 
2020 (Baden et al., 2020; Polack et al., 2020). However, under- 
standing vaccine-induced immunity is of high importance given 
the goal to achieve immunity for most people through vaccina- 
tion, rather than as a consequence of infection. 

The receptor binding domain (RBD) of the SARS-CoV-2 spike 
is an important target for serological and B cell studies because it 
directly interacts with the cellular receptor angiotensin convert- 
ing enzyme 2 (ACE2) which mediates host cell entry (Letko 
et al., 2020; Wrapp et al., 2020). Antibodies binding to the RBD 
can potently block attachment of the virus to ACE2 and thereby 
neutralize the virus (Barnes et al., 2020). As a consequence, 
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RBD-based vaccines are in development in addition to full- 
length spike-based vaccines (Krammer, 2020). Analyses of the 
B cell responses to the spike generally focus on the RBD and 
on cells sorted with RBD baits introducing an inherent bias by 
omitting non-RBD targets (Cao et al., 2020; Gaebler et al., 
2020; Robbiani et al., 2020; Weisblum et al., 2020). This is also 
true for B cells and monoclonal antibodies (mAbs) isolated 
from vaccinated individuals (Wang et al., 2021). However, other 
epitopes within the spike protein, notably the N-terminal domain 
(NTD) but also S2, do harbor neutralizing epitopes (Chi et al., 
2020; Liu et al., 2020; McCallum et al., 2021; Song et al., 
2020). In fact, the NTD is heavily mutated in the three most prom- 
inent variants of concern (VOCs) (B.1.1.7, B.1.351, and P.1) (Da- 
vies et al., 2021; Faria et al., 2021; Tegally et al., 2020). Here, we 
studied the unbiased plasmablast response to SARS-CoV-2 
mRNA-based vaccination and report several new findings. First, 
we document that RBD and NTD co-dominate as B cell targets 
on the viral spike protein, highlighting the importance of the 
NTD. We also report the first vaccine-induced NTD mAbs. In 
addition, we show that the majority of mAbs isolated are non- 
neutralizing, which is reflective of the higher binding to neutrali- 
zation ratios found in serum after vaccination compared to natu- 
ral infection. Finally, data from plasmablasts suggest that, at 
least some of the vaccine-induced response is biased by pre-ex- 
isting immunity to human B-coronaviruses. 


RESULTS 


The polyclonal antibody response to mRNA vaccination 
exceeds titers seen in convalescent individuals but is 
characterized by a high ratio of non-neutralizing 
antibodies 

In late 2020, six adult participants of an ongoing observational 
study received mRNA-based SARS-CoV-2 vaccines (Table 
S1). Blood from these individuals (termed V1-V6) was collected 
at several time points including before vaccination (for 4/6), after 
the first vaccination and at several time points after the second 
vaccination. We examined their immune responses to recombi- 
nant spike protein and RBD in enzyme-linked immunosorbent 
assays (ELISA), in comparison to those of 30 COVID-19 survivors 
(Figures 1A and 1B, Table S1). The sera from convalescent indi- 
viduals were selected based on their anti-spike titers and group- 
ed into three groups (low +: n = 8; moderate ++: n = 11; and 
high +++: n = 11, based on the antibody titer measured in the 
Mount Sinai’s CLIA laboratory (Wajnberg et al., 2020), taken 
111-273 days after symptom onset), in order to facilitate identi- 
fying different features that may track with the strength of the 
antibody response. Five out of six vaccinees produced anti- 
spike and anti-RBD responses that were, at the peak, markedly 
higher than responses observed even in the high titer convales- 
cent group while one vaccinee (V4) produced titers comparable 
to the high titer group. Notably, the antibody response peaked 
1 week after the second vaccine dose, followed by a decline in 
titers over the following weeks as expected from an antibody 
response to vaccination. Interestingly, anti-RBD antibody titers 
seemed to decline faster than anti-spike antibody titers, which 
appeared to be more stable over time. We also measured 
neutralizing antibody titers using authentic SARS-CoV-2 and 
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found a similar trend with all vaccinees displaying high titers, 
even though V4 responded with delayed kinetics (Figure 1C). 
Importantly, although at the peak response, the vaccine group 
mounted neutralization titers that fell in the upper range for the 
high convalescent group, they did not exceed that group mark- 
edly. This finding prompted us to calculate the proportions of 
spike binding to neutralizing antibodies. For the convalescent 
group, we found that individuals with lower titers had a higher 
proportion of binding to neutralizing antibodies than high-re- 
sponding convalescent individuals (Figure 1D). When deter- 
mined at the time of peak response, the vaccinees had the high- 
est proportion of binding to neutralizing antibody titers, 
indicating an immune response more focused on non-neutral- 
izing antibodies or an induction of less potent neutralizing anti- 
bodies in general (or both). These proportions remained stable 
over time with the ratio of binding to neutralizing antibodies in 
vaccinated individuals being significantly higher than those 
observed for any of the three convalescent groups (p = 0.0004, 
0.0002 and 0.0041 for the three groups respectively; 
Figure S1). We also investigated the spike binding to RBD bind- 
ing ratio and found no difference to convalescent individuals 
except a general trend toward proportionally less RBD binding 
over time in the vaccinees (Figure S1). 


mRNA vaccination induces a modest but measurable 
immune response to seasonal B-coronavirus spike 
proteins 

It has been reported that SARS-CoV-2 infection induces an orig- 
inal antigenic sin-type immune response against human corona- 
viruses (hCoVs) to which the majority of the human population 
has pre-existing immunity (Aydillo et al., 2021; Song et al., 
2020). Here, we explored whether this phenomenon is also 
induced by SARS-CoV-2 mRNA vaccination. Antibody titers in 
four vaccinees against spike protein from a-coronaviruses 
229E and NL63 were detectable at the pre-vaccination time 
point but did not increase substantially after vaccination (Figures 
1E and 1F; for V5 and V6 no pre-vaccination serum was avail- 
able). However, titers against the spike proteins of B-coronavi- 
ruses OC43 and HKU1 increased substantially in these four vac- 
cinees after vaccination (Figures 1G and 1H). Thus, vaccination 
with mRNA SARS-CoV-2 spike also boosts immune responses 
against seasonal B-coronavirus spike proteins in a manner remi- 
niscent of that reported for natural infection with SARS-CoV-2. 


The plasmablast response to SARS-CoV-2 mRNA 
vaccination targets both the RBD and the NTD 

In order to characterize the B cell response to vaccination in an 
unbiased manner, plasmablasts were single cell sorted from 
blood specimens obtained from three individuals (V3, V5, and 
V6) 6 days after the booster immunization (Figure S2). All mAbs 
were generated from single-cell sorted plasmablasts and probed 
for binding to recombinant SARS-CoV-2 spike protein. A total of 
21 (40 mAbs were screened, with 28 being clonally unique, Table 
S2) spike-reactive mAbs were isolated from V3, 6 (82 screened, 
20 unique) from V5, and 15 (84 screened, 24 unique) from V6 
(Figure 2A). Using recombinant spike, RBD, NTD, and S2 pro- 
teins, we mapped the domains to which these mAbs bind. Inter- 
estingly, only a minority of these antibodies recognized RBD 
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Figure 1. Antibody responses in individuals vaccinated with mRNA-based SARS-CoV-2 vaccines 

(A-C) Antibody responses of convalescent individuals and vaccinees to full-length spike protein (A) and RBD (B) as measured by ELISA and neutralizing activity of 
the sera of the same individuals in a microneutralization assay against authentic SARS-CoV-2 (C). Convalescent individuals were grouped based on their initial 
antibody response (measured in a CLIA laboratory) to spike protein into +, ++, and +++. 

(D) Ratios between binding and neutralizing antibody levels in vaccinees and convalescent individuals. Higher ratios indicate a bias toward non-neutralizing 


antibodies. 


(E-H) show antibody responses against -coronavirus 229E and NL63 and -coronavirus OC43 and HKU1 spike proteins over time. Bars represent the geometric 


mean, error bars represent the 95% confidence intervals. 


(24% for V3, 47% for V6, and no RBD binders were identified for 
V5) (Figures 2B and 2E). A substantial number of the isolated 
mAbs bound to NTD including 14% for V3, 33% for V5, and 
33% for V6 (Figures 2C and 2E). These data indicate that RBD 
and NTD are co-dominant in the context of the mRNA-induced 
plasmablast response. The epitopes for the majority of the re- 
maining spike binding mAbs, 52% for V3, 50% for V5, and 
20% for V6, mapped to S2 (Figures 2D and 2E). Only three 
mAbs were not accounted for in terms of binding target (two 
for V3 and one for V5; Figure 2E). 
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The majority of isolated mAbs from SARS-CoV2 
vaccinees are non-neutralizing 

All antibodies were tested for neutralizing activity against the 
USA-WA1/2020 strain of SARS-CoV-2. Only a minority of the 
binding antibodies, even those targeting the RBD, showed 
neutralizing activity (Figures 2F and 2G). For V3, only one (an 
RBD binder) out of 21 mAbs (5%) displayed neutralizing activity 
(Figure 2G). For V5, a single NTD antibody neutralized authentic 
SARS-CoV-2 (17%) (Figure 2G). The highest frequency of 
neutralizing antibodies was found in V6 (83%) with one RBD 
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Figure 2. Characterization of mAbs derived from vaccine plasmablasts 
(A-D) Binding of plasmablasts derived from three vaccinees (V3, V5, and V6) against full-length spike (A), RBD (B), NTD (C), and S2 (D). 


(E) The percentages of the respective antibodies per subject. 


(F and G) Neutralizing activity of the mAbs against authentic SARS-CoV-2 (F) and the proportion of neutralizing antibodies per subject is shown in (G). 

(H and |) Reactivity of mAbs to spike protein of human B-coronaviruses OC43 and HKU1. 

MBC, minimal binding concentration. All experiments except data shown in (H) and (I) were performed in duplicates and the mean of the duplicates is shown with 
standard deviation. For (H) and (I), a representative dataset from a singlet ELISA run is shown. 


neutralizer and four NTD neutralizers (Figure 2G). Interestingly, 
the highest neutralizing potency was found for mAb PVI.V5-6, 
an NTD binder followed by PVI.V6-4, an RBD binder. 

We also tested all mAbs for reactivity to the spike proteins of 
the four hCoVs: 229E, NL63, HKU1, and OC43. No antibody 
binding to the spike proteins of a-coronaviruses 229E and 
NL63 was found but we identified five mAbs (including three 
from V3, one from V5, and one from V6) that bound, to varying 
degrees, to the spike of OC43, which, like SARS-CoV-2, is a B- 
coronavirus (Figure 2H). Three mAbs showed strong binding 
(PVI.V3-8, PVI.V3-12, and PVI.V6-1), while PVI.V3-17 showed 
an intermediate binding phenotype and PVI-V5-1 bound very 
weakly. Three of these mAbs also showed binding to the spike 
of HKU1, another B-coronavirus. Of these, PVI.V6-1 showed 


only very weak binding while PVI.3-8 and PVI.3-12 had low min- 
imal binding concentrations (MBCs) indicating higher affinity 
(Figure 2l). 


The spike-reactive plasmablast response is dominated 
by IgG1* cells and is comprised of a mixture of cells with 
low and high levels of somatic hypermutation (SHM) 
Single-cell RNA sequencing (ScRNA-seq) was performed on 
bulk sorted plasmablasts from the three vaccinees (V3, V5, V6) 
to comprehensively examine the transcriptional profile, isotype 
distribution, and somatic hypermutation (SHM) of vaccine- 
induced plasmablasts. We analyzed 4,584, 3,523, and 4,461 sin- 
gle cells from subjects V3, V5, and V6, respectively. We first 
verified the identity of sequenced cells as plasmablasts through 
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Figure 3. Characterization of bulk sorted plasmablasts via single-cell RNA sequencing 


C) UMAP overlay of BCR sequences with confirmed spike binding activity. 


ee ee _ m 


A) Uniform manifold approximation and projection (UMAP) of scRNA-seq from bulk plasmablast with recovered BCR sequences (purple) or unrecovered (gray). 
B) UMAP overlay of percent of cellular population expressing MZB1, PRDM1, and XPB1. Hexbin equals 80 individual cells. 


D) Proportional composition of heavy chains genes in the spike binding sequences broken down by sample. 


(E) Comparison of nucleotide-level mutation frequency in immunoglobulin heavy chain variable (IGHV) genes between plasmablasts clonally related to spike- 
binding mAbs from SARS-CoV-2 vaccinees, plasmablasts sorted from PBMCs 1 week after seasonal influenza vaccination and found in vaccine-responding B 
cell clones, and naive B cells found in blood of an influenza vaccinee (left); and between plasmablasts from SARS-CoV-2 vaccinees found to be clonally related to 
spike-binding mAbs that were, respectively, cross-reactive and non-cross-reactive to human f-coronaviruses spike proteins (right). p values were generated 


using a two-sided Kruskal-Wallis test with Dunn’s post-test (left) or a Mann-Whitney U test (right). 


the combined expression of B cell receptors (BCRs) (Figure 3A) 
and that of the canonical transcription as well as other factors 
essential for plasma cell differentiation, such as PRDM1, XBP7, 
and MZB1 (Figure 3B). To identify vaccine-responding B cell 
clones among the analyzed plasmablasts, we used scRNA-seq 
to also analyze gene expression and V(D)J libraries from the 
sorted plasmablasts and clonally matched the BCR sequences 
to those from which spike-specific mAbs had been made. Using 
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this method, we recovered 332, 7, and 1,384 BCR sequences 
from the scRNA-seg data that are clonally related to the spike- 
binding mAbs derived from subjects V3, V5, and V6, respectively 
(Figure 3C). It is important to note here that we were not able to 
recover clonally related sequences for all of the mAbs that we 
cloned and expressed from each of the three vaccinees. 

We next examined the isotype and IgG subclass distribution 
among the recovered sequences. IgG1 was by far the most 
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Figure 4. Mapping of the amino acid substitutions and deletions onto the structure of the SARS-CoV-2 spike glycoprotein 


(A) Mutations of the three major variants of concern B.1.1.7, B.1.315, and P.1. 


(B) These mutations mapped onto the structure of the spike glycoprotein (model generated by superposition of PDB: 6MOJ and 7C2L) (Chi et al., 2020; Lan et al., 
2020). One RBD in the up conformation (red) is bound with ACE2 receptor (pink). The NTD is colored blue and the various amino acid substitutions are shown as 
yellow spheres. One spike protomer is shown in bold colors while the other two are colored white. 

(C) Competition between ACE2 and neutralizing RBD targeting mAbs PVI.V3-9 and PVI.V6-4 for binding to RBD. 

(D) BLI-measured binding affinities of the RBD mutants to ACE2, as well as the calculated fold change compared to wild type, are shown in the table on the right. 


dominant isotype in the three vaccinees (Figure 3D). Finally, we 
assessed the level of somatic hypermutation (SHM) among the 
mAbs-related sequences from the three subjects. We used 
the SHM levels observed in human naive B cells and seasonal 
influenza virus vaccination-induced plasmablasts that were pre- 
viously published for comparison (Turner et al., 2020). Spike- 
reactive plasmablasts from V3 and V6 but not V5 had accumu- 
lated SHM at levels that are significantly greater than those 
observed with naive B cells (Figure 3E, left). Strikingly, the 
SHM level among V6 plasmablasts was equivalent to those 
observed after seasonal influenza virus vaccination (Figure 3E, 
left). We reasoned that the high level of SHM among spike-reac- 
tive plasmablasts may be derived from those targeting 
conserved epitopes that are shared with human £-coronavi- 
ruses. Indeed, we found that the SHM level among clones that 
are related to cross-reactive mAbs was significantly higher 
than their non-cross-reactive counterparts (Figure 3E, right). 


Competition of RBD binding neutralizing mAbs with 
ACE2 and affinity of variant RBDs for human ACE2 

Two mAbs were identified as neutralizing and binding to RBD. We 
wanted, therefore, to test whether they competed with ACE2 for 
RBD binding. Concentration-dependent competition was indeed 
observed for both mAbs demonstrating that inhibition of ACE2 


binding is the mechanism of action of the two mAbs (Figure 4). 
Since we prepared RBD proteins of viral variants of concern for 
analysis of antibody binding (see below), we also wanted to 
assess the affinity of each variant RBD for human ACE2. Using 
biolayer interferometry (BLI), we measured association and disso- 
ciation rates of the N501Y RBD mutant (B.1.1.7 carries that muta- 
tion as its sole RBD mutation), Y453F, as found in mink isolates 
(Larsen et al., 2021), N439K, which is found in some European 
clades (Thomson et al., 2021), a combination of Y453F and 
N439K, E484K (part of B.1.351 and P.1) as well as for the 
B.1.351 and the P.1 RBDs for a recombinant version of human 
ACE2 (Figures 4A, 4B, and 4D). Almost all of the single and double 
mutations in RBD tested increased affinity to human ACE2. Spe- 
cifically, N501Y and Y453F combined with N439K increased affin- 
ity for human ACE2 by 5-fold (Figures 4D and S3). In contrast, 
E484K on its own decreased affinity by 4-fold. Of note, the 
B.1.851 RBD affinity for ACE2 was comparable to that of the 
wild-type RBD. These data were confirmed using an ELISA-based 
method which showed the same trends (Figure S4). 


Binding profiles of polyclonal serum and mAbs to RBDs 
carrying mutations found in viral variants of concern 
Next, we assessed binding of sera from vaccinated individuals, 
COVID-19 survivors, and mAbs derived from plasmablasts to 
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Figure 5. Binding and neutralization of SARS-CoV-2 variants 
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(A-C) Binding of serum samples from convalescent individuals (A), vaccinees (B), and vaccine-derived mAbs (C) to a panel of RBD mutants is shown. The red line 
in (A) indicates the average reduction. Dotted lines in (A) and (B) indicate 100%, the line with smaller dots in (C) indicated reactivity of the anti-his coating control. 
For vaccinees, late samples (V1 = d89, V2 = d102, V3 = d47, V4 = d48, V5 = 49, and V6 = 48) were assayed. 

(D) The spike mutations of virus isolate PV14252 mapped on a structure of the SARS-CoV-2 spike protein with ACE2 (model generated by superposition of PDB: 


6MOJ and 7C2L) (Chi et al., 2020; Lan et al., 2020). 


(E and F) The inhibitory effect of vaccine serum and vaccine-derived neutralizing antibodies on both wild-type SARS-CoV-2 and PV14252. 
(G) Neutralizing activity of the plasmablast-derived neutralizing antibodies against wild-type, B.1.1.7, and B.1.351 virus isolates. Of note, these comparative 
assays were always performed side by side but sets were run by different operators and on a different Vero cell clone as the neutralization assays shown in 


Figure 2. 


variant RBDs. Our panel of RBDs includes published mAb 
escape mutants, RBD mutants detected by the Mount Sinai Hos- 
pital’s Pathogen Surveillance Program in patients seeking care 
at the Mount Sinai Health System in NYC, as well as mutations 
found in viral variants of interest and variants of concern (Baum 
et al., 2020; Greaney et al., 2021b; Larsen et al., 2021; Thomson 
et al., 2021; Weisblum et al., 2020). Serum from convalescent in- 
dividuals showed strong fluctuations depending on the viral 
variant (Figure 5A). In general, single mutants E406Q, E484k, 
and F490K exerted the biggest impact on binding. However, 
complete loss of binding was rare and 2- to 4-fold reduction in 
binding was more common. Interestingly, almost all sera bound 
better to N501Y RBD (B.1.1.7) than to wild-type RBD (average 
129% compared to wild-type). Conversely, the B.1.351 RBD 
Caused, on average, a 39% reduction in binding. The impact 
was slightly lower for the P.1 RBD (average 70% binding 
compared to wild-type). For sera from the six vaccinated individ- 
uals, however, the highest reduction seen was only 2-fold for 
E406Q, N440K, E484K, and F490K (Figure 5B). Of note, the vac- 
cinees’ later samples (V1 = d89, V2 = d102, V3 = d47, V4 = d48, 
V5 = 49, and V6 = 48) were assayed to allow for some affinity 
maturation. The highest reduction observed for E484K, F484A, 
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B.1.351, and P.1 were also approximately 2-fold but this did 
not apply to all six vaccinees. Some vaccinees maintained bind- 
ing levels against these RBDs at levels comparable to wild- 
type RBD. 

RBD binding mAbs were also tested for binding to the same 
variants. In general, mAbs maintained binding levels within 2- 
fold of the binding seen with the wild-type RBD with some ex- 
ceptions. In fact, for most mAbs, no impact on binding was 
observed (Figure 5C) with the exception of PVI.V3-9, which lost 
binding to the RBD carrying F486A. Although there was a nega- 
tive impact on binding of several mAbs to the B.1.351 variant, 
binding was almost unaffected by the mutations in the P.1 
variant RBD. Only one mAb, PVI.V6-4, showed a drop in binding 
to P.1. 


Escape of an NTD and E484K mutant virus from 
polyclonal post-vaccination serum is negligible but NTD 
mutations significantly impact the neutralizing activity 
of NTD binding mAbs 

Through the Mount Sinai Hospital’s Pathogen Surveillance Pro- 
gram, we had access to the SARS-CoV-2 isolate PV14252 
(Clade 20C, Pango lineage B.1) that featured two mutations 
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(W64R, L141Y) and one deletion (A142-145) in the NTD as well as 
the E484K mutation in the RBD (Figure 5D). To determine the 
susceptibility of this virus variant to neutralization by post-vacci- 
nation serum, we performed microneutralization assays. Wild- 
type SARS-CoV-2 and PV14252 were tested in parallel to ensure 
that the assay setup for both viruses allowed comparison. We 
found a relatively minor impact when testing polyclonal sera 
from vaccinees for neutralizing activity (Figure 5E). The activity 
of sera from V2, V5, and V6 slightly increased while the activity 
for V1, V3, and V4 decreased. Next, we tested the seven neutral- 
izing mAbs that we isolated from plasmablasts. Consistent with 
their binding profiles in the variant RBD ELISA, the two RBD 
mAbs neutralized both viruses with comparable efficiency 
(Figure 5F). In fact, the activity of PVI.V3-9 increased slightly 
(Figure 5F). In stark contrast, all five anti-NTD antibodies 
completely lost neutralizing activity against PV14252 due to mu- 
tations present in the NTD of this viral isolate. 


B.1.1.7 and B.1.351 partially escape from plasmablast- 
derived neutralizing antibodies 

We also tested the neutralizing activity of the two RBD and the 
five NTD antibodies against the variants of concern B.1.1.7 
and B.1.351 (Figure 5G). Both variants contain deletions as 
well as mutations in the NTD. In addition, B.1.1.7 carries the 
N501Y RBD mutation and B.1.351 carries N417K, E484K, and 
N501Y mutations in the RBD (Figures 4A and 4B). The two 
RBD binding antibodies lost no (PVI.V6-4) or little (PVI.V3-9) 
neutralizing activity against B.1.1.7. However, PVI.V3-9 lost all 
activity against B.1.351 and the remaining neutralizing activity 
of PVI.V6-4 was low (but measurable). All but one (PVI.V6-11) 
NTD mAbs lost neutralizing activity against B.1.1.7 and all of 
them lost neutralizing activity against B.1.351 once more high- 
lighting the importance of changes in the NTD for the antibody 
activity. 


DISCUSSION 


Our knowledge of B cell responses to SARS-CoV-2 mRNA vacci- 
nation remains incomplete. We urgently need information about 
the nature of polyclonal vaccine-induced responses as well as 
unbiased, in-depth analyses of plasmablast responses. Our 
data provide important new insights into these responses in 
comparison with immune responses to natural infection. Indeed, 
SARS-CoV-2 infection results in a very heterogeneous antibody 
response to the spike protein in terms of antibody quantity. In 
contrast, MRNA vaccination appears to induce a high antibody 
response of relatively homogeneous titers. However, we also 
found that vaccinees generate more non-neutralizing antibodies 
than COVID-19 survivors resulting in a lower ratio of neutralizing 
to binding antibodies. These data were already apparent in the 
early phase clinical trials but remained unrecognized at the 
time (Walsh et al., 2020). Interestingly, low-titer convalescent 
serum had the highest relative amount of neutralizing antibodies, 
whereas the proportion of binding antibodies was increased in 
sera with higher measured antibody titers. The majority of plas- 
mablasts sampled after vaccination do, in fact, produce non- 
neutralizing antibodies. Two recent studies have performed a 
similarly unbiased plasmablast analysis for individuals naturally 
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infected with SARS-CoV-2 (Cho et al., 2021; Huang et al., 
2021). Of course, the antibody response after SARS-CoV-2 
infection is not only targeting the spike protein but several other 
proteins expressed by the virus. When accounting for spike 
binding only, these studies report proportions of 44% and 25% 
neutralizing antibodies (Cho et al., 2021; Huang et al., 2021). 
While plasmablast analysis is in general not quantitative (e.g., 
one clone per clonotype is selected, etc.), our analysis of post- 
vaccination plasmablasts found a lower number of neutralizing 
antibodies (17%). 

Future studies are needed to reveal the role of non-neutralizing 
antibodies in SARS-CoV-2 immune protection. Indeed, antibody 
functions other than neutralization have been shown to correlate 
with protection (Bartsch et al., 2021; Gorman et al., 2021; Scha- 
fer et al., 2021). The importance of absolute antibody titers and 
not ratios is underscored by the fact that post-vaccination 
neutralization titers were equal to or exceeded the titers found 
in the high responder convalescent group. 

Of the four seasonal CoVs that are widely circulating in hu- 
mans, B-coronaviruses OC43 and HKU1 have higher homology 
to SARS-CoV-2 spike. Vaccinated individuals mounted a 
response to spike proteins from OC43 and HKU1 but not to «-co- 
ronaviruses 229E and NL63. This phenomenon resembles the 
immune imprinting described in influenza virus immunology 
and has already been shown for natural infection with SARS- 
CoV-2 where a “backboost” to B-coronaviruses was also found 
(Aydillo et al., 2021; Song et al., 2020). A few of the mAbs isolated 
in our study had, indeed, such a cross-reactive phenotype. It re- 
mains unclear whether these antibodies, which target mostly S2 
epitopes, contribute to protection against SARS-CoV-2, OC43, 
or HKU1 infection. However, the cross-reactive epitopes of 
mAbs that do bind SARS-CoV-2, HKU1, and OC43 spikes could 
provide the basis for future pan-B-coronavirus vaccines. While it 
is likely the case that the B cells producing these mAbs come 
from recall responses and were initially induced by human B-co- 
ronaviruses (which is supported by serology and of course the 
extensive SHM that the mAbs show), they could hypothetically 
also be de novo induced antibodies. While this is probably 
not the case, we cannot exclude this possibility with our current 
data. 

Another interesting point we noted is the co-dominance of 
RBD and NTD. Previous analyses of B cell responses to SARS- 
CoV-2 mRNA vaccination focused on cells baited by labeled 
RBD (Wang et al., 2021). We, in contrast, took an unbiased 
approach to sort and clone plasmablasts in an antigen-agnostic 
manner. We found similar levels of NTD and RBD binders with 
many mAbs binding to epitopes outside the RBD and the NTD. 
In one vaccinee not a single RBD binding mAbs was isolated 
with the caveats that the overall number of mAbs derived from 
that individual were low and their polyclonal serum antibody re- 
sponses included RBD recognition. These data suggest that the 
NTD, which also harbors neutralizing epitopes, is—at least—as 
important as the RBD and warrants as much attention. In fact, 
five out of seven neutralizing antibodies isolated in this study 
bound to the NTD and only two targeted the RBD. Recent 
studies analyzing the plasmablast response after natural infec- 
tion have found a similar co-dominance of RBD and NTD (Cho 
et al., 2021; Huang et al., 2021) with one study reporting 59 
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mAbs targeting the RBD, 64 targeting the NTD, and 46 binding 
outside of RBD and NTD and the second study finding 10 RBD 
mAbs, 13 non-RBD S1 binding mAbs (strongly suggesting NTD 
binding), and 9 mAbs targeting S2. Interestingly, and in contrast 
to our findings, a recent deep mutational scanning paper with 
sera from mMRNA-1273 vaccinees found a very strong RBD- 
focused response (Greaney et al., 2021a). Further characteriza- 
tion of the mAbs obtained in our study showed a complete loss of 
neutralization against an authentic, replication-competent 
variant virus that harbored extensive changes in the NTD. All 
NTD mAbs also lost neutralizing activity against B.1.351 and all 
but one lost activity against B.1.1.7. These observations may 
explain why a reduction in neutralization against the viral variant 
of concern B.1.1.7 is seen in some studies despite the fact the 
N501Y substitution in the RBD of this variant does not signifi- 
cantly impact binding and neutralizing activity (Emary et al., 
2021). The key role of NTD as target for antibodies has recently 
also been shown using memory B cell-derived mAbs (McCallum 
et al., 2021). 

In addition, we assessed the impact of different RBD muta- 
tions on affinity toward human ACE2. Interestingly, N501Y 
increased the affinity by 5-fold. This increase in receptor binding 
affinity may contribute to the higher infectivity of B.1.1.7, which 
carries this mutation in its RBD. In contrast, introduction of 
E484K reduced the affinity by 4-fold which may explain why virus 
variants carrying only the E484K mutation have rarely spread 
efficiently, although viruses carrying E484K have been detected 
since the fall of 2020 in a handful of patients receiving care at the 
Mount Sinai Health System and have also been reported in 
immunocompromised patients (Choi et al., 2020). It is tempting 
to speculate that the N501Y mutation enables the acquisition 
of E484K without a fitness loss. In fact, the B.1.851 RBD, which 
carries N501Y and E484K (as well as N417K) showed binding to 
hACE2 that was similar to wild-type RBD. Recently, B.1.1.7 
variant strains carrying E484K, in addition to N501Y, have been 
isolated in the UK (PHE, 2021), providing evidence for the hy- 
pothesis that N501Y enables acquisition of mutations in the 
RBD that may be detrimental to receptor binding. However, 
recent expansion of B.1.526, a lineage also featuring E484K 
but without N501Y in New York City, suggests that this fitness 
loss may be overcome by other, yet uncharacterized, changes 
in the virus as well (Annavajhala et al., 2021; Lasek-Nesselquist 
et al., 2021). Interestingly, binding of convalescent sera to the 
N501Y RBD was also increased, suggesting that changes that 
increase affinity for the receptor may also increase affinity of a 
set of antibodies that may mimic the receptor. 

We also noted that the two neutralizing antibodies against the 
RBD showed some reduced binding to a mutant RBD carrying 
the E484K mutation while having similar or even increased 
neutralizing potency against a variant virus carrying the E484K 
mutation as the only change in its RBD. The reduced affinity of 
the E484K variant RBD for hACE2 could render the virus more 
susceptible to RBD binding mAbs. Thus, an antibody binding 
to the RBD may just be more effective in interfering with a low-af- 
finity as compared to a high-affinity RBD-hACE2 interaction. 
Increased affinity as an escape mechanism for viruses has 
been described in the past (Hensley et al., 2009; O’Donnell 
et al., 2012) and the converse mechanism could be at play here. 
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Whether or not the current vaccines will provide effective pro- 
tection against circulating and emerging viral variants of concern 
is an important question which has gathered a lot of attention in 
early 2021. Our data indicate that reduction in binding to the 
E484K and B.1.351 variant RBDs was minor (often only 2-fold) 
compared to reported reduction in neutralization (which ranges 
from 6- to 8-fold to complete loss of neutralization) (Cele et al., 
2021; Wibmer et al., 2021; Wu et al., 2021). Although not tested 
here, it is likely that the reduction in binding to full-length spike is 
even lower, given the many epitopes on the spike other than NTD 
and RBD. The maintenance of binding to a large degree 
observed in this study suggests that viral variants will have a mi- 
nor impact on serological assays which are currently in wide use 
for medical, scientific, and public health reasons. Binding, non- 
neutralizing antibodies have also been shown to have a protec- 
tive effect in many viral infections (Asthagiri Arunkumar et al., 
2019; DiLillo et al., 2014; Saphire et al., 2018) and may be a factor 
in the substantial residual protection seen in the Johnson & John- 
son and Novavax vaccine trials against B.1.351 in South Africa 
(Shinde et al., 2021). Production of non-neutralizing antibodies 
may also play a role in protection by mRNA vaccines after the 
first dose, as it is substantial and occurs during a time when 
neutralizing antibody titers are either very low or absent (Baden 
et al., 2020; Dagan et al., 2021; Polack et al., 2020). Finally, 
although some antibodies may lose neutralizing activity due to 
reduced affinity, they do still bind. Furthermore, B cells with 
these specificities potentially could undergo affinity maturation 
after exposure to a variant virus or a variant spike-containing 
vaccine, leading to high-affinity antibodies to variant viruses of 
concern. 

In summary, we demonstrate that the antibody responses to 
SARS-CoV-2 mRNA vaccination comprise a large proportion 
of non-neutralizing antibodies and are co-dominated by NTD 
and RBD antibodies. The NTD portion of the spike represents, 
thus, an important vaccine target. Since all viral variants of 
concern are heavily mutated in this region, these observations 
warrant further attention to optimize SARS-CoV-2 vaccines. 
Finally, broadly cross-reactive mAbs to B-coronavirus spike pro- 
teins are induced after vaccination and suggest a potential 
development path for a pan- B-coronavirus vaccine. 


Limitations of the study 

While our study characterizes the antibody response after 
SARS-CoV-2 mRNA vaccination in detail, it has several limita- 
tions. The first limitation is the small number of study partici- 
pants, which makes this study a qualitative rather than a quanti- 
tative study. Another limitation is the lack of plasmablast analysis 
of SARS-CoV-2-infected individuals. We have compared our 
data with published data from plasmablast analysis after 
SARS-CoV-2 infection but a side-by-side comparison would 
have been more accurate. We have also not included a longitu- 
dinal analysis of the convalescent sera in the study but feel that 
providing a wide range of time points and titer levels offsets this 
limitation to a certain degree. In addition, while the clones of 
crossreactive plasmablasts are likely derived from the memory 
compartment and have likely been initially induced by seasonal 
coronavirus infections, we lack the ultimate proof for that since 
we did not analyze pre-vaccination memory B cells and B cell 
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receptor sequences. Finally, the burning question of whether the 
abundant non-neutralizing antibodies do have a protective effect 
in vivo will need to be elucidated by follow-up studies. 
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IgD-PerCP-Cy5.5 (clone IA6-2) 
CD19-PE (clone HIB19) 
CD38-PE-Cy7 (clone HIT2) 
CD3-Alexa 700 (clone HIT3a) 
Streptavidin (HRP) 

Biotin Anti-6X His tag® antibody 
Anti-mouse IgG HRP 

Mouse anti-histidine antibody 
Anti-human IgG HRP 

Bacterial and virus strains 
SARS-CoV-2 (isolate USA-WA1/2020) 
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SARS-CoV-2 hCoV-19/South Africa/ 
KRISP-K005325/2020 


SARS-CoV-2 Isolate hCoV-19/England/ 
204820464/2020 


Biological samples 
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Blood samples 

Chemicals, peptides, and recombinant proteins 
Zombie Aqua 
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RNase Inhibitor, Murine 

Maxima H Reverse Transcriptase 
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Chromium Next GEM Single Cell 5’ Kit v2 
Library Construction Kit 

Chromium Next GEM Chip K Single Cell Kit 
Chromium Single Cell Human BCR 
Amplification Kit 

Dual Index Kit TT Set A 

SPRiselect Reagent Kit 

High Sensitivity D5000 ScreenTape 

High Sensitivity D5000 Reagents 

Bovine Serum Albumin 

SARS-CoV-2 RBD 


SARs-CoV-2 variant RBDs: E406Q, N417V, 
N439K, N440K, Y453F, E484K, F486A, 
N487R, F490K, Q493R, B.1.1.7, N439K/ 
Y453F, B.1.351 and P.1 

SARS-CoV-2 Spike (2P) 

OC43 spike 


SOURCE 


BioLegend 
BioLegend 
BioLegend 
BioLegend 
BioLegend 
BioLegend 
Abcam 

Abcam 
Rockland antibodies and assays 
Takara 
Millipore Sigma 


BEI 
PVI, Mount Sinai 
BEI 


BEI 


PVI, Mount Sinai 
PVI, Mount Sinai 


BioLegend 

Clontech 

New England BioLabs 
Thermo 

Thermo 

10X Genomics 

10X Genomics 

10X Genomics 

10X Genomics 


10X Genomics 

Beckman Coulter 
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Agilent 
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Krammer lab, Mount Sinai 
Krammer lab, Mount Sinai 


Krammer lab, Mount Sinai 
Krammer lab, Mount Sinai 
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Cat#348208 
Cat#302254 
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NR-52281 
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NR-54009 


NR-54000 


V1-V6, 30 convalescents 
V3, V5, V6 


Cat#423101 
Cat#635013 
Cat# M0314L 
Cat#EP0753 
Cat#15918014 
Cat#1000263 
Cat#1000190 
Cat#1000286 
Cat# 1000253 


Cat#1000215 
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Cat# 5067-5592 
Cat# 5067-5593 
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N/A 
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N/A 
N/A 
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HKU spike 

229E spike 

NL63 spike 

SARS-CoV-2 NTD 
SARS-CoV-2 S2 
Remdesivir 

Critical commercial assays 
Anti-human IgG Fc Biosensors 
Ni-NTA (NTA) Biosensors 
Deposited data 


Human B cell receptor repertoire data after 
seasonal influenza vaccination 


Experimental models: cell lines 
Expi293F 

Vero.E6 cells 

Expi293F 

Oligonucleotides 

oligo-dT23VN primer, random Hexamers 
15t PCR primers & nested PCR primers 
Gibson Cloning Primers 
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pCAGGS SARS-CoV-2 


pCAGGS SARS-CoV-2 SA variant spike 
pCAGGS SARS-CoV-2 UK variant spike 
pCAGGS SARS-CoV-2 RBD 
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OC43 2P spike plasmid 
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Seurat (v3.2.2) 
scRepertoire (v1.1.3) 


schex (v1.3.0) 
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Change-O v0.4.6 
SHazaM v1.0.2 
GraphPad Prism v9.0.2 
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Krammer lab, Mount Sinai 
Krammer lab, Mount Sinai 
Krammer lab, Mount Sinai 
Sino Biological 

Sino Biological 

NIH 


Sartorius Corporation 
Sartorius Corporation 


Turner et al., 2020 


GIBCO 
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Thermo Fisher 
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Integrated DNA Technologies 


https://www.ncbi.nim.nih.gov/pmc/ 
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This paper 
This paper 
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32398876/ 


This paper 
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articles/PMC5584442/ 
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articles/PMC5584442/ 


10x Genomics 
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Borcherding et al., 2020 


Saskia Freytag 
Ye et al., 2013 


Gupta et al., 2015 
Gupta et al., 2015 
GraphPad Software 
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N/A 


Cat# 18-5001 
Cat# 18-5101 
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51-01-13-03 
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N/A 


httos://support. 10xgenomics.com/ 
https://satijalab.org/seurat/ 


httos://github.com/ncborcherding/ 
scRepertoire 


https://github.com/SaskiaFreytag/schex 


https://ftp.ncbi.nih.gov/blast/executables/ 
igblast/release/1.14.0/ 


https://changeo.readthedocs.io/ 
httos://shazam.readthedocs.io/ 
httos://www.graphpad.com 
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REAGENT or RESOURCE SOURCE IDENTIFIER 

BLItz Pro 1.3.1.3 Forté Bio https://www.fortebio.com/blitz.html 

Microsoft Excel Microsoft https://www.microsoft.com/en-ww/ 
microsoft-365/excel 

Other 

International ImMunoGenetTics Information Giudicelli et al., 2005 http://www.imgt.org/vquest/refsegh. 

System (IMGT) human immunoglobulin html#refdir 


germline references, release 201931-4 


RESOURCE AVAILABILITY 


Lead contact 
Requests for information or reagents should be directed to Florian Krammer (florian.krammer@mssm.edu). 


Materials availability 

Plasmids for SARS-CoV-2 antigens have been deposited at BEI Resources and can also be requested from the authors. Plasmids for 
human coronavirus spikes can be requested from NIH’s Vaccine Research Center. MAbs and plasmids for mAb expression can be 
obtained from the authors upon reasonable request. Variant viruses can be sourced from BEI Resources. 


Data and code availability 
The published article contains all datasets analyzed during the study except for BCR sequencing data which can be requested from 
Ali H. Ellebedy upon reasonable request. 


EXPERIMENTAL MODEL AND SUBJECT DETAILS 


Human subjects and specimen collection 
The study protocols for the collection of clinical specimens from individuals with and without SARS-CoV-2 infection by the Person- 
alized Virology Initiative were reviewed and approved by the Mount Sinai Hospital Institutional Review Board (IRB-16-16772; IRB-16- 
00791; IRB-20-03374). All participants provided written informed consent prior to collection of soecimen and clinical information. All 
specimens were coded prior to processing and analysis. An overview of the characteristics of the vaccinees as well as the study par- 
ticipants with and without COVID-19 is provided in Table S1. The vaccinees received two doses of the Pfizer mRNA vaccine. 
Whole blood was collected via phlebotomy in serum separator tubes (SST) or ethylenediaminetetraacetic acid (EDTA) tubes. 
Serum was collected after centrifugation as per manufacturers’ instructions. Peripheral blood mononuclear cell (PBMCs) isolation 
was performed by density gradient centrifugation using SepMate tubes (Stemcell) according to manufacturers’ instructions. PBMCs 
were cryo-preserved and stored in liquid nitrogen until analysis. 


Recombinant proteins 

All recombinant proteins were produced using Expi293F cells (Life Technologies). Receptor binding domain (RBD) and spike protein 
of SARS-CoV-2 (GenBank: MN908947.3) was cloned into a mammalian expression vector, DCAGGS as described earlier (Amanat 
et al., 2020b; Stad|bauer et al., 2020). RBD mutants were generated in the pCAGGS RBD construct by changing single residues using 
mutagenesis primers. All proteins were purified after transient transfections with each respective plasmid. Six-hundred million Ex- 
pi293F cells were transfected using the ExpiFectamine 293 Transfection Kit and purified DNA. Supernatants were collected on 
day four post transfection, centrifuged at 4,000 g for 20 min and finally, the supernatant was filtered using a 0.22 um filter. Ni-NTA 
agarose (QIAGEN) was used to purify the protein via gravity flow and proteins were eluted as previously described (Amanat et al., 
2020b; Stadibauer et al., 2020). The buffer was exchanged using Amicon centrifugal units (EMD Millipore) and all recombinant pro- 
teins were finally re-suspended in phosphate buffered saline (PBS). Proteins were also run on a sodium dodecyl sulfate (SDS) poly- 
acrylamide gels (5%-20% gradient; Bio-Rad) to check for purity (Amanat et al., 2018; Margine et al., 2013). Plasmids to express re- 
combinant spike proteins of 229E, HKU1, NL63 and OC43 were generously provided by Dr. Barney Graham (Pallesen et al., 2017). 
NTD and S2 proteins were purchased from SinoBiologics. 


METHOD DETAILS 
ELISA 


Ninety-six well plates (Immulon 4 HBX; Thermo Scientific) were coated overnight at 4°C with recombinant proteins at a concentration 
of 2 ug/mL in PBS (GIBCO; Life Technologies) and 50 uls/well. The next day, the coating solution was discarded. One hundred uls per 
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well of 3% non-fat milk prepared in PBS (Life Technologies) containing 0.01% Tween-20 (TPBS; Fisher Scientific) was added to the 
plates to block the plates for 1 h at room temperature (RT). All serum dilutions were prepared in 1% non-fat milk prepared in TPBS. All 
serum samples were diluted 3-fold starting at a dilution of 1:50. After the blocking step, serum dilutions were added to the respective 
plates for two h at RT. Next, plates were washed thrice with 250 uls/well of TPBS to remove any residual primary antibody. Secondary 
antibody solution was prepared in 1% non-fat milk in TPBS as well and 100 uls/well was added to the plates for 1 h at RT. For human 
samples, anti-human IgG conjugated to horseradish peroxidase (HRP) was used at a dilution of 1:3000 (Millipore Sigma; catalog 
#A0293). For mouse samples, anti-mouse IgG conjugated to HRP was used at the same dilution (Rockland antibodies and assays; 
catalog #610-4302). Specifically, a mouse anti-histidine antibody (Takara; catalog #631212) was used as a positive control to detect 
proteins with a hexa-histidine tag. Once the secondary incubation was done, plates were again washed thrice with 250uls/well of 
TPBS. Developing solution was made in 0.05M phosphate-citrate buffer at pH 5 using o-phenylenediamine dihydrochloride tablets 
(Sigma-Aldrich; OPD) at a final concentration of 0.04 mg/mL. One hundred uls/well of developing solution was added to each plate for 
exactly 10 min after which the reaction was halted with addition of 50 uls/well of 3M hydrochloric acid (HCI). Plates were read at an 
optical density of 490 nanometers using a Synergy 4 (BioTek) plate reader. Eight wells on each plate received no primary antibody 
(blank wells) and the optical density in those wells was used to assess background. Area under the curve was calculated by deduct- 
ing the average of blank values plus 3 times standard deviation of the blank values. All data was analyzed in Graphpad Prism 7. This 
protocol has been described in detail earlier (Bailey et al., 2019; Wohlbold et al., 2015a). 

Purified monoclonal antibodies were used at a concentration of 30 ug/mL and then subsequently diluted 3-fold. Purified mono- 
clonal antibodies were only incubated on the coated plates for an h. The remaining part of the protocol was the same as above (Ama- 
nat et al., 2020a; Wohlbold et al., 2015b). 


Bio-layer Interferometry Binding Experiments 

Bio-layer Interferometry (BLI) experiments were performed using the BLItz system (forteBlO, Pall Corporation). Recombinant human 
Fc fusion ACE2 (SinoBiological) was immobilized on an anti-human IgG Fc biosensor, and RBDs were then applied to obtain binding 
affinities. Single-hit concentrations were tested at 5.8 uM for binding. All measurements were repeated in subsequent independent 
experiments. Kp values were obtained through local fit of the curves by applying a 1:1 binding isotherm model using vendor-supplied 
software. All experiments were performed in PBS pH 7.4 and at room temperature. 


hACE2 competition interferometry experiments 

Interferometry experiments were performed using a BLItz instrument (fortéBlO, Sartorius). Polyhistidine-tagged Fabs were immobi- 
lized on Ni-NTA biosensors at 10 ng/mL and SARS-CoV-2 RBD was supplied as analyte at 5uM alone or pre-mixed with hACE2-Fc at 
different concentrations. Maximal signal at association (Rmax) was used to plot the concentration-dependent competition with 
hACE2. All experiments were performed in PBS at pH 7.4 and at room temperature. 


RBD-hACE2 ELISA 

25ng of hACE2-Fc fusion protein expressed in HEK293 cells were adhered to high-capacity binding, 96 well-plates (Corning) over- 
night in PBS. Plates were blocked with 5% BSA in PBS containing Tween-20 (PBS-T) for 1 h at room temperature (RT). Blocking so- 
lution was discarded and 5-fold dilutions of 6xHis-tagged RBDs in PBS were added to wells and incubated for 1 h at RT. Plates were 
then washed three times with PBS-T. Anti-polyhistidine IgG-Biotin (Abcam) in PBS-T was added to each and incubated for 1 h at RT. 
Plates were then washed three times with PBS-T. Streptavidin-HRP (horseradish peroxidase) (Abcam) in PBS-T was added to each 
and incubated for 1 h at RT. Plates were then washed three times with PBS-T Plates were developed using 1-Step Ultra TMB 
(3,3’,5,5’-tetramethylbenzidine) substrate (ThermoFisher), stopped with sulfuric acid and immediately read using a plate reader at 
450nm. Data were plotted using Prism 9 (GraphPad Software) and affinities determined by applying a nonlinear regression model. 


Viruses and cells 

Vero.E6 cells (ATCC #CRL-1586) cells were maintained in culture using Dulbecco’s Modified Eagles Medium (DMEM, GIBCO) which 
was supplemented with 10% fetal bovine serum (FBS, Corning) and antibiotics solution containing 10,000 units/mL of penicillin and 
10,000 ug/mL of streptomycin (Pen Strep, GIBCO)(10). Wild type SARS-CoV-2 (isolate USA-WA1/2020), hCoV-19/South Africa/ 
KRISP-K005325/2020 (B.1.351, BEI Resources NR-54009) and hCoV-19/England/204820464/2020 (B.1.1.7, BEI Resources NR- 
54000) were grown in cells for 3 days, the supernatant was clarified by centrifugation at 4,000 g for 5 min and aliquots were frozen 
at —80°C for long term use. The viruses were subjected to deep sequencing to ensure that no mutations had taken place in culture. 
The polybasic cleavage site changed to WRAR in the B.1.351 variant virus during cultivation in cell culture (as known for this virus at 
BEI Resources) and no other unexpected mutations occurred. A primary virus isolate, PV14252, bearing mutations and deletions in 
the spike was obtained by incubating 200 uls of viral transport media from the nasopharyngeal swab with Vero.E6 cells. The 
sequence of the passage 2 viral isolate was identical to the sequence obtained directly from the clinical specimen. Sequencing 
was performed on the Illumina platform as described previously (Gonzalez-Reiche et al., 2020). The replication competent viruses 
were used to test serum from study participants and antibodies for neutralization activity. 
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Neutralization assay 

Twenty-thousand cells in 100 uls per well were seeded on sterile 96-well cell culture plates one day prior to the neutralization assay. In 
general, cells were used at 90% confluency to perform the assay. All serum samples were heat-inactivated to eliminate any comple- 
ment activity. Serial dilutions of serum samples were made in 1X minimal essential medium (MEM; Life Technologies) starting at a 
dilution of 1:20. All work with authentic SARS-CoV-2 (isolate USA-WA1/2020 and PV14252) was done in a biosafety level 3 (BSL3) 
laboratory following institutional biosafety guidelines and has been described in much greater detail earlier (Amanat et al., 2020b; 
Amanat et al., 2020c). Six hundred median cell culture infectious doses (TCIDs 9s) of authentic virus (USA-WA1/2020 and 
PV14252) was added to each serum dilution and the virus-serum mixture was incubated together for 1 h inside the biosafety cabinet. 
Media from the cells was removed and 120 uls of the virus-serum mixture was added onto the cells for 1 h at 37°C. After one h, the 
viruS-serum mixture was removed and 100 uls of each corresponding dilution was added to every well. In addition, 100uls of 1X MEM 
was also added to every well. Cells were incubated for 48 h at 37°C after which the media was removed and 150 uls of 10% form- 
aldehyde (Polysciences) was added to inactivate the virus. For assay control, remdesivir was used against both the wild type virus as 
well as the patient isolate. After 24 h, cells were permeabilized and stained using an anti-nucleoprotein antibody 1C7 as discussed in 
detail earlier (Amanat et al., 2020b; Sun et al., 2020). 


Cell sorting and flow cytometry 

Staining for sorting was performed using cryo-preserved PBMCs in 2% FBS and 2 mM ethylenediaminetetraacetic acid (EDTA) in 
PBS (P2). Cells were stained for 30 min on ice with CD20-Pacific Blue (2H7, 1:400), Zombie Aqua, CD71-FITC (CY1G4, 1:200), 
IgD-PerCP-Cy5.5 (IA6-2, 1:200), CD19-PE (HIB19, 1:200), CD38-PE-Cy7 (HIT2, 1:200), and CD3-Alexa 700 (HIT3a, 1:200), all Bio- 
Legend. Cells were washed twice, and single plasmablasts (live singlet CD19* CD3° IgD'° CD38* CD20" CD71") were sorted using 
a FACSAria Il into 96-well plates containing 2 uL Lysis Buffer (Clontech) supplemented with 1 U/uL RNase inhibitor (NEB) and imme- 
diately frozen on dry ice, or bulk sorted into PBS supplemented with 0.05% BSA and processed for single cell RNaseq. 


Monoclonal antibody (mAb) generation 

Antibodies were cloned as described previously (Wrammert et al., 2011). Briefly, VH, Vk, and VA genes were amplified by reverse 
transcription-PCR and nested PCR reactions from singly sorted plasmablasts using primer combinations specific for IgG, IgM/A, 
Igx, and Igà from previously described primer sets (Smith et al., 2009) and then sequenced. To generate recombinant antibodies, 
restriction sites were incorporated via PCR with primers to the corresponding heavy and light chain V and J genes. The amplified 
VH, Vk, and VA genes were cloned into IgG1 and lgx expression vectors, respectively, as described previously (Nachbagauer 
et al., 2018; Wrammert et al., 2008). Heavy and light chain plasmids were co-transfected into Expi293F cells (GIBCO) for expression, 
and antibody was purified with protein A agarose (Invitrogen). 


Single-cell RNA-seg library preparation and sequencing 

Bulk-sorted plasmablasts were processed using the following 10x Genomics kits: Chromium Next GEM Single Cell 5’ Kit v2 (PN- 
1000263); Library Construction Kit (PN-1000190); Chromium Next GEM Chip K Single Cell Kit (PN-1000286); Chromium Single 
Cell Human BCR Amplification Kit (PN- 1000253), and Dual Index Kit TT Set A (PN-1000215). The cDNAs were prepared after 
GEM generation and barcoding, followed by GEM RT reaction and bead cleanup steps. Purified cDNA was amplified for 10-14 cycles 
before cleaning with SPRiselect beads. Then, samples were evaluated on a 4200 TapeStation (Agilent) to determine cDNA concen- 
tration. B cell receptor (BCR) target enrichments were performed on full-length cDNA. Gene expression and enriched BCR libraries 
were prepared as recommended by the Chromium Next GEM Single Cell 5’ Reagent Kits v2 (Dual Index) user guide, with appropriate 
modifications to the PCR cycles based on the calculated cDNA concentration. The cDNA libraries were sequenced on Novaseg S4 
(Illumina), targeting a median sequencing depth of 50,000 and 5,000 read pairs per cell for gene expression and BCR libraries, 
respectively. 


Single-cell RNA-seq analysis 

Single-cell RNA sequencing and BCR sequencing data was processed using Cell Ranger v5.0 and the GRCh38-2020 version of the 
human genome provided by the manufacturer. Total recovered cells by RNA sequencing were V3: 6,608, V5: 5,256, and V6: 6,325 
with a mean of 90.64% read mapped to the genome. Count matrices were processed in R (v4.0.2) using the Seurat (v3.2.2) R package 
(Stuart et al., 2019). Cells were filtered for percentage of mitochondrial genes less than 15% and number features less than 4,000. The 
three specimen sequencing runs were integrated using log-normalized count values and canonical correlation approach (Stuart et al., 
2019) with 2,000 variable features. The resulting single-cell object underwent principal component analysis and the top 30 principal 
components were used for uniform manifold approximation and projection and identifying neighbors. Clustering was performed us- 
ing aresolution of 0.6. The integrated RNA sequencing object included 12,568 cells with V3: 4,584, V5: 3,523, and V6: 4,461 cells. The 
filtered contig annotation outputs of Cell Ranger vdj were loaded into R and processed using the scRepertoire (v1.1.3) R package 
(Borcherding et al., 2020). Clonotypes were assigned using igraph (v1.2.6) network analysis of components generated from CDR3 
sequences greater than or equal to 0.85 normalized Levenshtein distance. Percent of cells expressing genes along the UMAP 
embedding was visualized using the schex (v1.3.0) R package. For mutation analysis, heavy chains of mAbs and single-cell BCRs 
first underwent V(D)J gene annotation using IgBLAST (v1.14.0) (Ye et al., 2013) with human reference (release 201931-4) from the 
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international ImMunoGenefTics information system (IMGT) (Giudicelli et al., 2005) and then parsing using Change-O (v0.4.6) (Gupta 
et al., 2015). Mutation frequency was calculated, as described in (Turner et al., 2020), using the “calcObservedMutations” function 
from SHazaM (v.1.0.2) (Gupta et al., 2015) and by counting the number of nucleotide mismatches from the germline sequence in the 
heavy chain variable segment leading up to the complementary-determining region 3 (CDR3), while excluding the first 18 positions 
that could be error-prone due to the primers used for generating the mAb sequences. 


QUANTIFICATION AND STATISTICAL ANALYSIS 


Structural figures were modeled and rendered in PyYMOL (The PyMOL Molecular Graphics System, Version 2.4 Schrodinger, LLC). 
Statistical analysis was performed in GraphPad Prism using a one-way ANOVA with correction for multiple Comparisons for 
Figure S1. For Figure 3, p values were generated using a two-sided Kruskal-Wallis test with Dunn’s post-test or a Mann-Whitney 
U test. Significance was defined as p < 0.05, p values are directly indicated in graphs. Number of subjects can be found in the results 
and methods section, definition of center, and dispersion and precision measures are described in the respective figure legends. 
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Figure S1. Comparison of binding to neutralizing titer ratios between naturally infected and vaccinated individuals (A) and full-length spike to 
RBD ratios (B), related to Figure 1 


Statistical analysis was performed in GraphPad Prism using a one-way ANOVA with correction for multiple comparisons, significance was defined as p < 0.05. 
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Figure S2. Gating strategy for sorting plasmablasts from total PBMCs isolated 1 week after second immunization, related to Figure 2 and 3 
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Figure S3. Representative Bio-Layer Interferometry binding isotherms from two independent experiments, related to Figure 4 
The raw data are shown in pink and the Langmuir 1:1 kinetics fit is shown in black. 
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Figure S4. Binding of SARS-CoV-2 variant RBDs to ACE2, related to Figure 4 


(A) ELISA curves of the RBD variants binding to human ACE2. Shown are the binding curves calculated with nonlinear regression to the arithmetic mean values 
from eight replicates + SEM. The calculated steady-state Kp values + SEM from end-point ELISA measurements and the fold-change in comparison to wild type 
RBD are reported in (B). 





